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Abstract. Modern control and diagnostic systems (CDS) usually determine only the technical 
condition (TC) at the current time, ie the CDS answers the question: a complex technical system 
(CTS) should be considered operational or not, and may provide little information on performance 
CTS even in the near future. Therefore, the existing scenarios of CDS operation do not provide for 
the assessment of the possibility of gradual failures, ie there is no forecasting of the technical condi- 
tion. 

The processes of parameter degradation and degradation prediction are stochastic processes, 
the “behavior” of which is influenced by a combination of external and internal factors, so the deg- 
radation process can be described as a function that depends on changes in the internal parameters 
of CTS. 

The hybrid method involves the following steps. The first is to determine the set of initial 
characteristics that characterize the CTS vehicle. The second is the establishment of precautionary 
tolerances of degradation values of the characteristics that characterize the pre-failure technical 
condition of the CTS. The third is to determine the rational composition of informative indicators, 
which maximally determine the "behavior" of the initial characteristics. The fourth — implementa- 
tion of multiparameter monitoring, fixation of values of the controlled characteristics, formation of 
an information array of values of characteristics. Fifth — the adoption of a general model of the 
process of changing the characteristics of the CTS. Sixth — the formation of a real model of the 
process of changing the characteristics of Y(t) on the basis of an information array of values of 
characteristics obtained by multi-parameter monitoring. Seventh — forecasting the time of possible 
occurrence of the pre-failure state of the CTS, which is carried out by extrapolating the obtained 
real model of the process of changing the characteristics of Y(t). It is proposed to use two types of 
models: for medium- and long-term forecasting - polynomial models, for short-term forecasting — 
a linear extrapolation model. 

At the final stage, forecast errors are determined for all types of models of degradation of pa- 
rameters and characteristics. Based on the results of the forecast verification, the models are adjust- 
ed. 

Keywords: technical condition, parameter degradation, multiparameter monitoring, pre- 
failure state, polynomial models, extrapolation model. 
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Anotania. CygacHi CHcCTeMH KOHTposIFO Ta WiarHoctyBaHHa (CKJ]) Bu3HayaroTb, 3a3BM4all, JIM- 
Ie TexHIYHuii ctaH (TC) y norouHuii MoMeHT Yacy, TOOTO CK] BizmoBiqae Ha MMTaHHA: CKIayHy Tex- 
HigHy cuctemy (CTC) ciliy BBaxKaTH Tpalle3qaTHOrO 4H Hi, 1 IPH WbOMy MOxKe HaaBaTH OOMAalIb BIO- 
MOCTei LOO MpatesqaTHocti CTC Hasitb y HaiiOMKGOMy MaliOyTHBOMy. Omke, HasBHi ClleHapii po- 
doTu CK] ne nepexOayaroTb OLHIOBaHHA MOXKIMBOCTI BAHHKHCHHA MOCTYMOBUX Bi{MOB, TOOTO Hemae 
IIPOrHO3yBaHHA TeXHI4HOTO CTaHy. 

IIpouecn yerpayauii napaMertpis Ta IporHo3yBaHHA Aerpayauii — We croxacTH4Hi Mpouecn, 
Ha “TlOBeIHKy” AKHX BIVIMBa€ CYKYIMHICTb 30BHIIMHIX Ta BHYTPiHIX (akTOpiB, TOMy Tpowec e- 
rpaflaiii MoxKHa OMMCaTH AK (PYHKUIO, WO 3aJ1eKUT Bi 3MiHM BHYTPIWIHIx Mapametpis CTC. 

TidpuqHHui MeToy MepexOayae Taki etanu. lepmummii — Bu3HadeHHaA HaOopy BUXIHUX Xapak- 
TepHCTHK, AKi XapakTepu3yroTh TC CTC. Jpyruii — BcTaHoBseHHa MomepeKYBaIbHUX JOIMYCKIB 
B3Ha4eHb Jlerpayaiii xapakTepHCTHK, AKi XapakKTepH3yIOTb epexBINMOBHUM TexHi4HHi craH CTC. 
Tpetiii — Bu3Ha4YeHHA pallioHaIbHoro cKayy iHPOpMaTHBHHX MOKa3HHKIB, AKl MAaKCHMAaJIbHO 
OOYMOBJIIOIOTS “TIOBeEIHKy” BUXiHHX XapakTepucTuK. UeTBepruli — 3fiicHeHHa OaraTomapamet- 
PHYHOFO MOHITOPHHTY, (Pikcallid 3HAaYCHb KOHTPOJIbOBAaHUX XapaKTepHCTHK, (opMyBaHHA indop- 
MalllMHOTO MacuBy 3HaveHb xapakTepucTuK. IT’ aruii — npHiiHATTA 3arabHOi Moses Mmpolecy 3MiHu 
xapaxtepuctuk CTC. Woctuii — dopmysanna peasbHoi Moses Mpolecy 3MiHH XapaKTepHCTHK 
Y(t) Ha 0a31 inbopMaliifHOro MacuBy 3Ha4eHb XapaKTepHCTHK, 3100yYTOTO WWIAXOM 3/liicHeHHA Oara- 
TOMapaMeTpH4HOrO MOHITOpHHTy. Cbhomuli — MmporHo3yBaHHA MOMCHTY 4acy MOXKIMBOI MOABH Te- 
peaBiqMoBHoro craHy CTC, ake 3iCHIOETBCA WIIAXOM eKCTpallouAii 3400yTo! peabHoi MOJeui 
lIpolecy 3MiHu xapaktepucTuK Y(t). [IpononyeTbca BHKOPHCTOByBAaTH Ba BUM Mojesei: WIA cepe- 
J{HBO- Ta JJOBFOCTPOKOBOrO TIPOrHO3yBaHHA — TOJHOMIAJIbHI MOJ{eII, Jd KOPOTKOCTPOKOBOrO TIpo- 
THO3yBaHHA — JNHIMHYy eKCTpanoAMitHy MOeJb. 

Ha 3aBeplasIbHOMY eTalli 30iMCHIO€TBCA BH3HAYCHHA MOMMJIOK TIPOrHo3y AIA ycix BUIB MO- 
fee erpafaii MapaMerpis Ta xapakTepucTHuK. 3a pe3yIbTaTaMU Bepudikalli mporHosy 3Ailic- 
HIOETBCA KOPHTyBaHHA MOJesIeH. 


Ks1r040Bi C1OBa: TeXHIWHH cTaH, Werpayauli Mapametpis, OaraTomapaMeTpH4Huii MOHITO- 
PHHT, MepexBiTMOBHUH CTaH, MOHOMIAJIbHI MOJeII, EKCTpanloaAiitHa MOesIb. 


Modern control and diagnostic (CDS) performance does not cover the full range of 
systems usually determine only the technical possible failure situations during their classi- 
condition (TC) of complex technical systems fication. Learning systems are nonlinear, they 
(CTS) at the current time. The existing scenar- are able to classify situations even in condi- 
ios of CDS operation do not provide for fore- tions of uncertainty, but they do not provide 
casting the technical condition. the necessary reliability. 

An analysis of recent research and pub- It is proposed to develop a hybrid meth- 
lications on the assessment of the current od of diagnosis, which is based on a combina- 
technical condition and its forecasting shows tion of approach, which allows to perform al- 
that almost all known studies are based on the so medium- and long-term forecasts, and the 
assumption of the stationary process of for- approach to perform, including operational 
mation of gradual failures, which is erroneous. forecasting. The hybrid method involves sev- 

Degradation processes and degradation en steps. 
prediction are stochastic processes. Their "be- At the final stage, forecast errors are de- 
havior" is influenced, first, by a combination termined for all types of models of degrada- 
of external and internal factors; secondly, the tion of parameters and characteristics. Ac- 
degradation process can be described as a cording to the results of the forecast verifica- 
function that depends on changes in the inter- tion (provided that the errors exceed the limit 
nal parameters of the CTS. As a result, the use values that are set in advance), the adjustment 
of traditional approaches, which require a of, first, the models obtained with the help of 
very large knowledge base, leads to high both MGUA and BMLE is carried out; sec- 


computational complexity. Insufficient CDS 


fe 
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ondly, monitoring parameters, in particular 
time quantization intervals. 


Statement of the problem in general 

Assessment of the state of complex 
technical systems (CTS) is usually carried out 
using a monitoring and diagnostic system 
(MDS). However, modern CDS usually de- 
termine only the technical condition (TC) at 
the current time, CDS answers the question: 
CTS should be considered operational or not, 
and may provide little information about the 
operability of CTS even in the near future. 
Therefore, the existing scenarios of CDS do 
not provide for the assessment of the possibil- 
ity of failures, which are called gradual, there 
is no forecasting of the technical condition [1] 
— [4], which allows to predict the failure of 
CTS and timely initiate appropriate preventive 
actions to maintain serviceability of CTS. 

Thus, the issues of correct formulation 
of the procedure for assessing and forecasting 
the performance of STS with the use of mul- 
tiparameter monitoring, the formation of a list 
of indicators by which it becomes possible to 
carry out current assessment, as well as accu- 
rate and reliable forecasting, should be con- 
sidered relevant. 


Analysis of recent research and pub- 
lications. An analysis of recent research and 
publications on the assessment of the current 
technical condition and its forecasting [1] — 
[6] shows that much attention has been paid to 
determining the current technical condition of 
technical facilities, but much less research has 
been carried out related to TC forecasting. 
Moreover, almost all known studies are based 
on the assumption of the stationary process of 
formation of gradual failures, which is erro- 
neous and, accordingly, a significant disad- 
vantage of such methodological approaches. 

Another area of research is the use of 
accumulated information about the values of 
current parameters that determine the TC, 
which becomes possible due to almost con- 
stant multi-parameter monitoring. However, 
as a rule, forecasting of the TC is either not 
carried out, reducing the study to the construc- 
tion of the current TS signature and its valida- 
tion, or obtain a forecast with large errors [7] 


— [9]. 
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In addition, it should be noted that the 
approach, which is based on the use of accu- 
mulated information about the values of cur- 
rent parameters, requires significant compu- 
ting power - especially the speed of infor- 
mation processing and the establishment of 
primitive causal relationships. This is due to 
the fact that the application to implement the 
model of degradation of the TC, for example, 
the method of group argumentation requires, 
first, a significant sample; secondly, such a 
polynomial model is not capable of sensitivi- 
ty, ie does not take into account the "rate" of 
changes in parameter values and, accordingly, 
does not allow regulation of the duration of 
observation intervals (according to, for exam- 
ple, Kotelnikov's theorem) and prediction. in 
determining and forecasting the TC, and its 
size. In addition, polynomial models, which 
are created by the method of group argumen- 
tation requires, do not have such a quality as 
emergence, ie the model is not able to take 
into account the sudden factors that affect the 
controlled indicator. As a result, it can be con- 
cluded that the use of method of group argu- 
mentation requires (and other similar meth- 
ods) may be appropriate, primarily to create a 
medium- and long-term forecast [3]. 

It is worth noting that attempts have al- 
ready been made to determine the technical 
condition of the CTS and predict it. An exam- 
ple is the work [10], which proposes to use 
the Berg method and Shi models based on 
Bayesian networks to solve problems of diag- 
nosing and forecasting of complex technical 
systems. But most researchers still believe 
that the most convenient device for solving 
problems of diagnosis and prediction are arti- 
ficial neural networks [11] — [15]. 

Thus, it can be reasonably concluded 
that the methodological approaches currently 
used to determine the TS are not sufficiently 
consistent with the solution of the problem of 
forecasting the technical condition. 


The purpose of the article 

The purpose of the article is to present a 
hybrid method of intelligent diagnosis and 
prediction of complex technical systems using 
multi-parameter monitoring, as well as correc- 
tion of monitoring parameters during observa- 
tions. 
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Presentation of the main research mate- 
rial. 

First of all, it should be determined that 
the processes of degradation and, accordingly, 
the prediction of degradation, which are the 
subject of this article, are stochastic processes, 
the "behavior" of which is first influenced by 
a combination of external and internal factors; 
secondly, the degradation process can be de- 
scribed as a function that depends on changes 
in the internal parameters of the CTS. As a 
result, the use of traditional approaches, which 
require a very large knowledge base, leads to 
high computational complexity. Insufficient 
CDS performance does not cover the full 
range of possible failure situations during 
their classification. Learning systems are non- 
linear, they are able to classify situations even 
in conditions of uncertainty, but they do not 
provide the necessary reliability [9]. Thus, 
given the fact that none of the above ap- 
proaches to the diagnosis of CTS can be used 
alone, automatic diagnosis of the TC based on 
observations obtained during monitoring re- 
quires the development of a specific hybrid 
method. It is proposed to develop a method of 
diagnosis, which is based on a combination of 
approach, which allows to perform also medi- 
um- and long-term forecasts, and the approach 
to perform, including operational forecasting. 

The hybrid method involves the follow- 
ing steps [16] — [18]: 

a) determination of a set of initial char- 
acteristics {Y1, Y2, ..., Yn}, which character- 
ize the TC of CTS with a predetermined prob- 
ability. Output characteristics can be parame- 
ters of output signals (for example, power, 
frequency, pulse duration, pulse amplitude, 
etc.) or some system parameters (amplitude 
and phase characteristics, bandwidth, sensitiv- 
ity, etc.); 

b) establishment of precautionary toler- 
ances of degradation values of characteristics 
{Y1, Y2, ..., Yn}, which characterize the pre- 
failure technical condition of CTS. Advance 
tolerance is understood as the value of the 
characteristics between the limit and_pre- 
failure levels. The output of the value of the 
characteristic for the limit level means failure, 
and the achievement of the pre-failure level 
— the need for maintenance or replacement of 
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elements of the CTS, ie inconsistencies of the 
first or second type; 

c) determination of the rational compo- 
sition of informative indicators that determine 
the maximum "behavior" of the initial charac- 
teristics: 

{Yi}=f{Xn}, 


where K=(x, Xq ye Xy — "input" pa- 


rameters; 

d) implementation of multiparameter 
monitoring, fixation of values of controlled 
characteristics, formation of an information 
array of values of characteristics, which 
should be presented in the form of a matrix: 


yi ij 
Tage) 2 Ss ane 
yl mn 
where: 
m — is the total number of controlled 
characteristics of the CTS; 
n — is the number of measurements of 


each characteristic on the monitoring inter- 
val T. 

In the resulting matrix, the number of 
columns increases as monitoring. Each col- 
umn represents a point in the Euclidean con- 
figuration space, which should be denoted by 
{i = 1, ..., n}, and each row is a sequence of {j 
= /, ..., m} measurements of the values of the 
controlled characteristics; 

d) adoption of a general model of the 
process of changing the characteristics of 
CTS — function Y(t)={Yi(t)}"i=1,, and it 
should be borne in mind that the peculiarity of 
the process of changing the characteristics of 
Y (t) is its pronounced non-stationary 
throughout the monitoring interval, due to the 
presence in Y (t) of an irreversible component 
that characterizes the aging or wear process; 

e) formation of a real model of the pro- 
cess of changing the characteristics of Y (t) on 
the basis of an information array of values of 
characteristics obtained by multi-parameter 
monitoring. As a rule, the most adequate 
models are statistical, based on time series; 

f) forecasting the time of possible occur- 
rence of the pre-failure state of the CTS, 
which is carried out by extrapolating the ob- 
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tained real model of the process of changing 
the characteristics of Y (t). 

The procedure for forecasting the 
technical condition is to form a model of 
some a posteriori (conditional on the ob- 
served implementation of Y (t)) random pro- 
cess based on monitoring data and a priori 
information and further evaluate its charac- 
teristics. The formed model of a posteriori 
process characterizes individual properties 
of CTS. The result of the forecast, at the 
same time, gives an assessment of the con- 
dition, which is valid for the observed CTS. 

It is proposed to use two types of mod- 
els: 

—for medium- and long-term forecast- 
ing — polynomial models, which are ob- 
tained using the well-known MSUA [19]. 
But it should be noted that as basic models 
are widely used not only polynomials, but 
also nonlinear, probabilistic functions or 
clustering [20]; 

—for short-term forecasting — a linear 
extrapolation model, which is obtained using 
the method of multidimensional linear extrap- 
olation (MDLE), the use of which allows to 
perform the task of forecasting the technical 
condition of the CTS. 

The essence of the MDLE method is 
quite traditional finding the values of indica- 
tors of a system at a point outside the observa- 
tion segment, their values within this segment. 
The problem of creating an approximation 
model (A-Model) is a problem of spatial in- 
terpolation and extrapolation, which is de- 
scribed by such terms. 

Let there be a finite set of experiments- 
points in the space of experiments in the g-th 
period of time T, where some indicators of the 
system are defined. The task is to estimate the 
value of the vector of indicators of the system 
in a new experiment that is not contained in 
the specified set at the q-th time t under the 
condition: t, ¢ T. 

In general, the problem is formulated as 
follows. It is necessary to find the transfor- 
mation operator F from the equation: 

Y =F (X), 


where X= {X7, X2, ..., Xs, t} 1s the vector of 
input parameters that determines the coordi- 
nates of the experience point in the subspace 
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of parameters and on the time axis, which is 
the Euclidean space of dimension S; 

S — the number of parameters on which 
the characteristics of the system depend; 

Y= {Y1, Y2, ..., Ym, t} is a characteristic 
vector that determines the coordinates of the 
experience point in the characteristics sub- 
space, which is a Euclidean space of dimen- 
sion m; 

m is the number of system characteris- 
tics; 

F is a transformation operator that maps 
a subspace of parameters to a subspace of 
characteristics. 

The problem is to recover the vector F 
by some set of vectors {X}, for which the 
values of the vectors {Y} are known. That is, 
when observing a real system, a training sam- 
ple is obtained — a number of K points for 
which the vector Y is assigned to each vector 
X at time t. 

The use of the MDLE method assumes 
that the process of changing the values of pa- 
rameters and characteristics is piecewise line- 
ar throughout the direct use of CTS. It follows 
from this assumption: 

Subspaces of parameters and character- 
istics are linear. Thus, the equations for the 
elements of these subspaces can be written as 
follows: 


k=l 
{X} = Xi + S Ai (Xi+1 — X1), 


i=l 


k-1 
(Yp=¥it+>) wi (Vins -Y1), 


i=l 


where X,, i=1k — vectors in the subspace 


of parameters, each of which corresponds to a 
vector in the subspace of characteristics 


Y,i =1,k on the basis of k measurements; 


Xi, Hi, 1=1,S1,k — coefficients of the 
corresponding vector linear subspaces. 

2. The mapping of the subspace of pa- 
rameters to the subspace of characteristics is 
linear — Ai = pi. 

Thus, the task of forecasting is reduced 
to the following stages: sampling; calculation 
of the predicted value of the characteristic Y 
*: determination of forecast errors. 

Based on the results of observing the 
change in the parameters X and _ the 
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characteristics Y, measured at time intervals Tt, 
a sample of k (k2S+J/) implementations of the 
S-dimensional vector is formed. The interval t 
of fixing the values of X; and Y; is chosen so 
that if the forecast states that the value of the 
characteristic Y* exceeds the tolerance a, B, 
the advance t would allow to take 
organizational and technical measures to 
prevent failure. 

From k measuring points k'(k' <S + 1) 
of implementations which satisfy functions 


O(X, 4 ) which is a criterion of proximity is 
chosen: 
—_ pee —= 
D(X,, X= (X, — Xi’, O(X,, X) > min, 
i=l 
where Pa is the value of the i-th parameter 


that lies outside the observation segment. 

By k 'measurement points, the function 
F is restored. 

The next step is to calculate the linearity 
coefficients 4. To do this, a system of 
equations: 


s s Ss s 
Ati +A, > highs $2 tA YD hile a >i —U, 
i=l i=l i=l i=l 


S S s S 
Ay iilta +A, hi tc) alia = as —U, 
i=l i=l i=l i=l 


Vex, =H 1=1,8: 


hy=Xj+1i-—X1i « 
The system is written as an equation using a 
matrix form: 
H'’HA=H'U, (1) 
Nig iigassl igs) A, 
erie Mais en Miguse ee As 
Taig lies Avg 
U, 
u=|°? 
Us 


The coefficient A is from equation (1): 
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A= HU ("= Hy". 


The predicted value of the characteristic Y * 
is determined from the expression: 


kil 
Y=Yi+ > 4 (%a-K), 
i=l 


where wi= Ai, i =1,k according to the 
previously accepted assumption. 

At the final stage, forecast errors are 
determined for all types of models of 
degradation of parameters and 
characteristics. According to the results of 
forecast verification (provided that the 
errors exceed the limit values (@, @), which 
are set in advance), the adjustment is 
obtained, first, the models obtained using 
both MSUA and = BMLE; © secondly, 
monitoring parameters, in particular time 
quantization intervals At; (Fig. 1). 


Ato 


Ati At; Atiss t 


Fig. 1. The procedure for determining the discreteness 
of fixing parameter values 


To solve this problem, we hypothesize 
that the process of changing the values of the 
observed parameter €(t) is stationary. In this 
case, the value of the discreteness of the time 
of fixing the values of the process of changing 
the parameters of the CTS is determined 
based on the expected model of the process of 
changing the values of the parameter &(t) 
(usually the time series is selected as a model 
of such a process) [21]. To form a model, ac- 
cording to the results of monitoring the pro- 
cess of changes in the values of the parameter 
€(t), an approximation function is constructed. 
Obtaining such a function in analytical form 
makes it possible to calculate it, after which, 
taking into account the value of the upper lim- 
it frequency of the spectrum of the approxi- 
mation function, according to Kotelnikov's 
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theorem, the value of the time discreteness of 
fixing the values of parameters At is calculat- 
ed. It should be emphasized that when the 
condition of stationarity (quasi-stationarity) of 
the process on the observation interval of the 
value At=const. 

However, the real processes of changing 
the values of the parameters of the CTS are 
often not stationary and in this case may re- 
quire correction of both the model and, ac- 
cordingly, the value of the discreteness of the 
fixation time. To do this, at each fixation 
point #; (tj; ¢ 7, T is the monitoring interval of 
the parameter €(t)), the values of the parame- 


ter € (t) are calculated =. If 
model of the process of changing the values of 
the controlled parameter is not adjusted, and, 
accordingly, on the observation interval the 
value At = const. If the condition is fulfilled 
after the next control step, which is performed 
at the point tj, then the spectrum of the func- 


dg 
dt 
ed again and, accordingly, the discreteness of 
the observation interval At; + 1 is calculated. 
Thus, the discreteness of the observation is 
adjusted. 

After adjusting the model in the next 
measurement, the stationary of the process of 
changing the values of the parameter € (t) is 
checked — the second derivative of the func- 
tion — the model of the process of changing 
the values of the parameters is calculated 
é"(t). If the condition is met ¢"(t)=0, the 
model is not adjusted € 0, if it is necessary to 
reform the model of the process of changing 
the values of € (t) [22]. 

A simplified block diagram of the pro- 
posed methodological approach to assessing 
and forecasting the technical condition of the 
CTS is presented in Fig. 2. 


“| =0, then the 
dt 


tion 


[| 0, its upper frequency is calculat- 


Conclusion 

It is proposed a hybrid method of 
intelligent diagnostics and forecasting of the 
condition of complex technical systems with 
the use of multiparameter monitoring, as well 
as correction of monitoring parameters during 
observations of the technical condition of 
complex technical systems (CTS). The 
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method also provides for the possibility of 
adjusting the type or model of functions that 
describe the processes of changing the values 
of parameters and characteristics of the 
technical condition. Adjustment is possible 
and appropriate in case of non-stationary 
process observed. 
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START 


The results of 
parameters 
control 


Calculation of the value of 
residual of the current val- 
ues of the parameters and 
the values of the parame- 
ters in the front step T; 


Collection and processing 
of measured values of pa- 


rameters with pace T : ; 
Creating an approximate 


Calculation of the value of model (A-Model) 


the residual parameters 


Normative values of the 
parameters and residual 


Calculation of the magni- 


Comparison of parameter tude of the residual param- 
values with their normative eA 


values 


The current 
value of the 

parameters in 
step Ti 


Comperison of the values 
of the predicted parameters 
with their normative values 


Determination of the state Forecasting parameters 
of working capacity values 


Exponential results of 
monitoring and prognostic 
parameters 


FINISH 


Fig. 2. Simplified block diagram of the algorithm for diagnosing and forecasting the technical condition 
of a complex technical system 
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